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ABSTRACT 

The formation of quasar black holes during the hydrodynamic collapse of protogalactic 

gas clouds is discussed. The dissipational collapse and long-term dynamical evolution of 

these systems is analysed using three-dimensional numerical simulations. The calculations 

0\ focus on the final collapse stages of the inner baryonic component and therefore ignore the 

j^ presence of dark matter. Two types of initial conditions are considered: uniformly rotating 



^ 



spherical clouds, and irrotational ellipsoidal clouds. In both cases the clouds are initially cold, 
^~~^ homogeneous, and not far from rotational support (r/|H^| ^ 0.1). Although the details of 

\Q the dynamical evolution depend sensitively on the initial conditions, the qualitative features 

Q of the final configurations do not. Most of the gas is found to fragment into small dense 



o 

0^ 



clumps, that eventually make up a spheroidal component resembling a galactic bulge. About 

5% of the initial mass remains in the form of a smooth disk of gas supported by rotation in 

f^ the gravitational potential well of the outer spheroid. If a central seed black hole of mass 

I ^ 10 M© forms, it can grow by steady accretion from the disk and reach a typical quasar 

O 



black hole mass ~ 10 M© in less than 5 X 10 yr. In the absence of a sufficiently massive 
seed, dynamical instabilities in a strongly self-gravitating inner region of the disk will inhibit 
steady accretion of gas and may prevent the immediate formation of a quasar. 
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1. INTRODUCTION 

The earliest objects observed to exist in the universe are quasars (Schneider, Schmidt, & 
Gunn 1991). Their bright emission is often associated with the accretion of gas onto massive 
(10 M0) black holes. Current cosmological simulations (Katz et al. 1993) predict the 
formation of potential quasar sites at redshifts as high as 2; = 8. The corresponding density 
peaks involve massive (~ 10 Mq), self-gravitating clouds of cold gas that have collapsed 
to a scale ^ 1 kpc, below the spatial resolution of the simulations. However, it is far from 
clear whether this cold gas can collapse further to form a massive black hole. In particular, 
the angular momentum of the gas will prevent radial infall beyond the centrifugal barrier, at 
a radius ~ 0.1-1 kpc. This is still larger than the Schwarzschild radius of the system by about 
6-8 orders of magnitude. The dynamics of the collapse is further complicated by small-scale 
fragmentation, and possibly also star formation and supernova explosions. In this paper, 
we study the three-dimensional hydrodynamics of the collapse and fragmentation of these 
rotating protogalactic gas clouds. We focus on hydrodynamics exclusively, concentrating 

on the possibility that massive black holes can form directly during the initial collapse of 
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protogalaxies at high redshifts . 

The existence of massive black holes is supported by various observational constraints on 
the compactness and efficiency of the energy source in active galactic nuclei. Recent observa- 
tions of unresolved gravitationally lensed images (indicating a source size ^ 1 pc from HST 
imaging of the quasar 2237+0305; see Rix, Schneider & Bahcall 1992), milliarcsecond jets 
(showing an unresolved core ^ 10 cm in some objects; see Baath et al. 1992), gravitational 
microlensing (indicating a continuum source of size ;^ 2 X 10 cm in 2237+0305; see Ranch 
& Blandford 1991), and rapid X-ray variability (Remillard et al. 1992), all favor the black 
hole paradigm over its starburst alternative (Terlevich et al. 1992 and references therein). 
The basic question is then: How did massive black holes form at high redshifts? The answer 
is obviously linked to the formation of structure in the early universe (Turner 1991). The 
coincidence between the peak in the quasar population at z ^ 2 and the expected epoch of 
galaxy formation in the cold dark matter (CDM) cosmology may be interpreted as evidence 
that a quasar activity reflects the initial stage in the formation of galactic bulges (Efstathiou 
& Rees 1988; Rees 1992; Haehnelt & Rees 1993). Fueling of quasar black holes may also have 



3 Alternatively, black holes could form as a result of the dynamical evolution of dense star clusters (Rees 
1984; Quinlan & Shapiro 1990). Even in this case, the final stages of the evolution will likely involve 
gas dynamics. 



been triggered by interactions between neighboring galactic systems (Hernquist 1989; Carl- 
berg 1990). However, there is presently no direct observational evidence indicating whether 
quasar black holes formed before or after galaxies were assembled (Loeb 1993a). Relatively 
young galaxies in the local universe do not show evidence for the level of activity found in 
high-redshift quasars. 

The central theoretical issues that need to be addressed with respect to the black hole 
formation process are as follows: 

1. Angular momentum barrier. Collapsing gas clouds acquire angular momentum 
through tidal interactions with neighboring perturbations (Hoyle 1949; Peebles 1969; Ef- 
stathiou & Jones 1979; White 1984; Hoffman 1986; Barnes & Efstathiou 1987; Ryden 1988; 
Quinn & Binney 1992). As the dark matter virializes, the gas cools and sinks to the center 
of the potential well until it becomes rotationally supported on scales of order ~ 0.1 — Ikpc. 
A black hole may form only if the specific angular momentum of the gas is reduced by 3-4 
orders of magnitude. The efficient loss of angular momentum can in principle be mediated 
by coUisional friction of free electrons or dust particles with the cosmic background radiation 
(Loeb 1993b; Umemura, Loeb & Turner 1993), but only for perturbations that collapse at 
very high redshifts [z ^ 100). Alternatively, quasars may be associated with rare systems 
that aquire low values of angular momentum from tidal torques during their cosmological 
collapse (Eisenstein & Loeb 1993). 

2. Fragmentation and star formation. As the gas cools and settles into a rotationally 
supported configuration the Jeans mass drops down to a stellar mass scale. Instead of 
being accreted by a central point mass the gas may then fragment into dense clumps or 
form stars. There is direct observational evidence that considerable fragmentation precedes 
quasar activity. Analysis of the broad emission lines in quasar spectra indicates a metallicity 
higher than solar (Hamann & Ferland 1992). An early phase of star formation and metal 
enrichment must therefore have occured in the quasar host systems. 

3. Formation of a single black-hole rather than a black hole cluster. If massive 
black holes resulted from the growth of many (e.g. stellar mass) seeds by accretion then 
black hole clusters were more likely to appear rather than single black holes. A black-hole 
cluster forming on the scale of a galactic bulge would extend much beyond the radius of 
the broad line region (~ lO^GM/c^ = O.lpc(M/lO^M0); cf. Peterson 1993) and would be 
in confiict with constraints on the compactness of the continuum source in quasars. The 
growth of massive black holes in bulges must therefore start from a limited number of seeds. 



For a primordial power spectrum of density fluctuations P{k) oc k'^ with n > —3, structure 
forms flrst on small scales and black holes could have formed directly near the Jeans mass 
(~ 10 M0). Abundant formation of black holes in many low mass systems (Gnedin & 
Ostriker 1992) that subsequently combine to make larger systems would result in black hole 
clusters. In contrast, the CDM power spectrum has a slope n Ri 3 at the low-mass end. The 
formation of black-hole clusters may then be avoided because all mass scales below 10 Mq 
collapse at about the same time. In addition, low-mass systems have shallower potential 
wells out of which the gas can easily be expelled by supernovae driven winds (Haehnelt & 
Rees 1993). 

4. What determines the minimum black-hole mass? The masses of central black holes 
in active galactic nuclei can be deduced observationally by a variety of methods: reverbera- 
tion mapping of the broad-line region (Peterson 1993), X-ray variability studies (Wandel & 
Mushotzky 1986), Eddington limit considerations (Haehnelt & Rees 1993), continuum flt of 
accretion disk models (Laor 1990), and photoionization models of the emission line intensi- 
ties (Netzer 1990; Wandel & Yahil 1985; Padovani, Burg & Edelson 1990). In all cases, the 
deduced black-hole masses are found to be ^ 10 Mq. Since the lowest mass objects are not 
necessarily observed to be close to their detection threshold, this lower bound on their mass 
needs to be explained. 

5. Fueling the central blak hole. Fueling of a quasar black hole requires high overdensities 
{Sp/p ^ 10 ) within a relatively massive host object (^ 10 Mq in gas) at high redshifts 
(Turner 1991). If a quasar radiates near the Eddington limit the radius of the fueling 
system should be ~ 0.6kpc(iv/10^^ erg • s~-^)-^'^(e/0.1)^'^ and its gas density must be ~ 
1.1 X 10 (e/0.1)~ Mq ■ kpc~ , where e is the conversion efficiency of the accreted mass into 
radiation and L is the quasar luminosity. Both of these characteristics correspond roughly 
to an evolved galactic bulge at high redshifts. It is therefore natural to link the activity in 
quasars to the collapse and evolution of galactic bulges. 

This work will address all of the above issues in the context of a simplified hydrodynamic 
model for the collapse of protogalactic gas clouds in the early universe. The discussion 
is complementary to previous studies that concentrated on the dynamics on larger scales 
(Efstathiou & Rees 1988; Katz 1992; Katz et al. 1993). At galactic bulge densities the 
cooling time of the optically thin gas is much shorter than the hydrodynamical time, and 
the gas temperature remains low Tg ^ 10 K. Under these conditions, pressure forces are 
negligible (except in shocks where they provide dissipation), and only rotation can provide 



support against gravity. Here we explore the hydrodynamics of protogalactic gas clouds in 
this regime using hybrid TV-body and smoothed particle hydrodynamics (SPH) simulations. 
The numerical method and the initial conditions for the simulations are described in §2. The 
results are presented in §3. In §4 we discuss the possible formation of a central seed black 
hole and its growth by accretion from a gaseous disk formed during the collapse. Finally, §5 
summarizes our main conclusions. 



2. NUMERICAL APPROACH AND BASIC ASSUMPTIONS 

2.1. Smoothed Particle Hydrodynamics 

The smoothed particle hydrodynamics (SPH) method has been used for the calculations 
presented here. SPH is a Lagrangian method that was introduced specifically to deal with as- 
trophysical problems involving self-gravitating fluids moving freely in three dimensions. The 
key idea is to calculate the pressure-gradient forces by kernel estimation, directly from the 
particle positions, rather than by flnite differencing on a grid, as is done in more traditional 
methods. This idea was introduced originally by Lucy (1977) and Gingold & Monaghan 
(1977), who applied it to the calculation of dynamical instabilities in rapidly rotating stars. 
Since then, a wide variety of astrophysical fluid dynamics problems have been tackled using 
SPH (see Monaghan 1992 for a recent review). In the past few years, these have included 
large-scale structure formation in cosmology (Katz, Hernquist, & Weinberg 1992) galaxy 
formation (Katz 1992), star formation (Monaghan & Lattanzio 1991), supernova explosions 
(Herant & Benz 1992), and stellar interactions (Rasio & Shapiro 1991, 1992). 

We use a modifled version of the SPH code developed originally by Rasio (1991) for the 
study of hydrodynamic stellar interactions (Rasio & Shapiro 1991, 1992). The implemen- 
tation of the SPH scheme is similar to that adopted by Hernquist & Katz (1989), but the 
gravitational fleld is calculated using a fast grid-based FFT solver. The neighbor searching 
in the code is also performed using a grid-based algorithm. Speciflcally, a multigrid hierar- 
chical version of the linked-list method usually adopted in P M particle codes (Hockney & 
Eastwood 1988) has been developed. The improved algorithm is extremely efficient, even 
for very nonuniform distributions of particles, provided that one is careful to fine-tune the 
ratio Lg/{hi)g of grid separation Lg to the average SPH smoothing length {hi)g for that 



grid. Other details about the implementation, as well as the results of a number of test-bed 
calculations were presented in the papers cited above. 

The most important modification made to the code for this study is the addition of a 
particle-mesh (PM) integrator that allows us to treat a coUisionless component in the system, 
in addition to the fluid component. The coUisionless particles are treated basically like SPH 
particles, except that only the gravitational force is retained in their equation of motion. A 
simple cloud-in-cell scheme (see, e.g., Hockney & Eastwood 1988) is used for both assigning 
density to, and interpolating forces from grid points. Thus the code can now be used for 
hybrid PM/SPH simulations. In the calculations of dissipational collapse described here, the 
PM method is used to model the system of compact high-density clumps of gas that form 
due to the gravitational instability. Speciflcally, whenever the density pi at the position of 
an SPH particle becomes larger than some critical value, the SPH particle is transformed 
instantaneously into a coUisionless particle. The gas represented by that SPH particle is 
assumed to continue its collapse to smaller and smaller scales, so that it effectively decouples 
from the rest of the fluid, except for the gravitational interaction. This procedure, although 
admittedly fairly crude, allows us to model the dynamical evolution of a collapsing system 
well past the point where the rapidly increasing dynamic range in fluid densities would 
make a purely hydrodynamic calculation grind to a halt (because the timestep must satisfy 
At ^ (G/Omax)~ )• The feedback effects of star formation on the large-scale hydrodynamics 
is ignored. Indeed, in the galaxy formation simulations of Katz (1992), which incorporated 
crudely the effects of star formation and energy production in supernovae, this feedback was 
found to have little effect on the numerical results. 

Another important modiflcation is that the SPH energy (or entropy) equation appro- 
priate for adiabatic conditions has been replaced by an isothermal equation of state. This 
is equivalent to setting the adiabatic exponent F = 1 and letting the "entropy variables" 
Ai = Pil p\ = Cg = k^Tg/ fi be constant for all particles (where pi is the pressure and p. the 
mean molecular weight). The assumption of a constant temperature is amply justifled in our 
model, since the gas is optically thin and cooling is so efficient that the temperature remains 
always ^ 10 K. The exact value of the temperature is unimportant, since pressure-gradient 
forces remain always negligible during the collapse. Shock dissipation is provided by the 
standard artificial viscosity, with the quadratic (von Neumann- Richtmyer) term dominating 
in this isothermal regime. 

All calculations shown here were done using AT % 10 particles (total, both SPH and 



coUisionless), and each SPH particle interacting with a nearly constant number of neighbors 
Nj^l ^ 64. The gravitational potential is calculated by FFT on a 256 grid. A typical run 
takes about 10 CPU hours on an IBM ES-9000 supercomputer. 

2.2. Initial Conditions 

For simplicity we consider an initially homogeneous cloud, represented by a random spa- 
tial distribution of SPH particles all having the same individual mass. We do not introduce a 
spectrum of small-scale density fluctuations, and avoid relaxing the initial distribution of SPH 
particles, since this would introduce strong correlations between particle positions on scales 
comparable to the particle smoothing lengths hi. We adopt units such that G = M = R = 1 
where M = Mg X 10 Mq is the total mass and R = R^ kpc is the initial size of the system. In 
these units the hydrodynamical time tdyn = (-R^/G-M)-^" Ri 10 jt R^' Mg~ ' . We assume 
an isothermal equation of state Pi = c^pi with c^ = 0.01 in our units, corresponding to a 
temperature Tg ^ b X 10 KMgi?^ [fi/mp), where rup is the proton mass. This is justifled 
since the cooling time of the gas tcool ^ 10 yr remains much smaller than the dynamical 
time t(jyn throughout the evolution. 

Angular momentum in the system is added in one of two ways. The flrst type of initial 
condition we consider is a rigidly rotating spherical cloud. The angular velocity is given by 
Q = 0.3 in our units, corresponding to a ratio of rotational kinetic energy to gravitational 
binding energy r/|H^| = ^0,"^ /{GM/R^) = 0.1 initially. The total angular momentum is 
then given by Jl{GM^Rfl'^ % 0.2. 

The second type of initial condition we consider is an irrotational (zero-vorticity) elliip- 
soidal cloud. To produce angular momentum without vorticity, we use the velocity fleld of a 
classical irrotational Riemann ellipsoid of type S (Chandrasekhar 1969, Chap. 7). Riemann 
ellipsoids are equilibrium fluid conflgurations supported by a combination of solid body ro- 
tation and internal fluid motions of uniform vorticity. In the type-S ellipsoids, the vorticity 
vector is everywhere parallel to the rotation axis. The irrotational ellipsoids are a particular 
case where the solid body rotation and internal motions combine to produce zero vorticity as 
seen in the inertial frame. Our initial condition is of course not in equilibrium, but the veloc- 
ity fleld is the same as in the equilibrium solution. In addition, we adjust the axis ratios of 
the ellipsoidal cloud so that the (uniform) density and total angular momentum are exactly 
the same as in the uniformly rotating sphere considered previously. The motivation here is 



that precisely this type of irrotational velocity field tends to be produced by external tidal 
torques acting on a mass of fiuid (see, e.g., Kochanek 1992). Vorticity cannot be produced 
by tidal torques in the absence of dissipation. Thus, an irrotational configuration is most 
likely if the angular momentum of the baryonic component in the protogalaxy was induced 
by external tidal torques before the collapse became dissipational. If, on the other hand, 
orbit crossing and dissipative interactions (e.g., with a dark matter halo; see Katz 1992) 
have already occurred, then the fiuid velocity field may have acquired vorticity and more 
organized fiuid motions resembling solid body rotation are possible. The two types of initial 
conditions considered in this paper are meant to represent these two different regimes. 



3. NUMERICAL RESULTS 

3.1. Collapse of a Uniformly Rotating Sphere 

The initial phase of the dynamical evolution (Fig. la, b) is characterized by the overall 
collapse of the system to a thin disk and the formation of dense filaments and clumps of 
cold gas. In the absense of an appropriate treatment, this process of continued fragmen- 
tation to smaller and smaller scales would make the numerical calculation grind to a halt 
well before a final equilibrium state is reached. This is because the hydrodynamic time 
inside clumps rapidly becomes orders of magnitude shorter than that of the system as a 
whole. As explained in §2, we solve this problem by converting gas particles into point-like 
objects whenever the local density at the particle's position is above a critical value /Omax- 
These point-like objects, representing high-density gas clumps at the limit of the numerical 
resolution, are treated subsequently as a coUisionless system interacting with the remaining 
gas though gravitational forces only. In the simulation of Fig. 1, we used pma,x/p = 10 , 
where p is the mean mass density in the system. This choice of /Omax is dictated primarily 
by practical considerations, since the total computation time of each simulation increases 
like /Omax- Calculations performed with different values of pma,x/ P in the range 10-10 show 
small differences in the details of the evolution, but remain all in good qualitative agreement. 

In just a few dynamical times, the dense gas clumps and filaments collide and merge, 
leading to fragmentation on smaller and smaller scales until a significant fraction, typically 
~95%, of the total mass is converted into the coUisionless component (Fig. 2). At interme- 
diate times (Fig. Ic, d), a tightly bound, compact group containing a small number (^ 10) 



of individual subsystems is observed. Further interactions between these subsystems lead to 
mergers, so that only two subsystems remain after a few crossing times (Fig. le), and these 
finally merge to produce the single structure obtained at the end of the calculation (Fig. If). 

This final, quasi-equilibrium configuration consists of a thin disk of cold gas embedded in 

4 
a virialized, coUisionless halo of a spheroidal shape . The gas is supported by rotation in 

the external halo potential. No measurable deviation from axisymmetry is left in this final 

configuration, although some evidence for the formation of a bar in the innermost region of 

the gaseous disk may be identified in Figure If (and is expected theoretically, cf. §4.1 below). 

We did not attempt to follow numerically the subsequent dynamical evolution of the gas in 

this inner region, since it is at the limit of our spatial resolution. 

3.2. Collapse of an Irrotational Ellipsoid 

The initial phase of the collapse of an irrotational ellipsoidal cloud is strikingly different 
(Fig. 3a, b). Here the system as a whole collapses to a spindle-like configuration at first. 
However, fragmentation on small scales and the formation of high-density gas clumps proceed 
as before. After several dynamical times, a small number of orbiting subsystems are again 
seen, although the overall mass distribution appears more centrally concentrated (Fig. 3c, 
d). The subsystems quickly merge into a single configuration (Fig. 3e,f): a virialized halo of 
dense gas clumps forms, while a small fraction ^ 4% of the mass remains in a thin gaseous 
disk. Thus, quite independent of the details of the initial conditions, it appears that the 
late stages of the dynamical evolution are always dominated by violent relaxation of the 
high-density clumps, and the collisions and mergers between a small number of subsystems. 
These complex dynamical processes tend to erase any memory of the initial conditions, and 
ultimately determine the principal characteristics of the final configuration. 



4 It is interesting to note that a similar structure is seen in the bulge of our own Galaxy, where radio 
observations show the presence of a cold disk containing ~ IO^Mq of molecular gas (i.e., a few percent 
of the total bulge mass) extending out to a distance ~ 300 pc from the Galactic center (Blitz et al. 
1993). 



4. DISCUSSION 



4.1. Disk Evolution 



It is tempting to envision that viscous accretion of the ~ 10 Mq of gas in the central 
disk onto a central seed object could lead to the formation of a supermassive black hole. 
This would occur on a viscous timescale, 

tvis ~ 5 X 10^ jt: mI^^ Rl^^a-\cs/10kms-^)-^, (1) 

where Cg is the sound speed in the disk, and a ^ 1 is a dimensionless measure of effective 
viscosity (Shakura & Sunyaev 1973). As the disk shrinks in size it may reach a point where 
its further contraction is limited by cooling rather than by angular momentum transport. In 
the process of losing angular momentum the disk radiates its gravitational binding energy 
with a luminosity, L ~ GM / Rty[s- This luminosity would approach the Eddington limit 
at a radius ~ O.lpc X Mg a ' (cs/lOkm • s~ ) ' . The subsequent evolution would still be 
limited by a relatively short cooling timescale, 

icool ~ ei^ < 4.4 X lOVr, (2) 

where e is the overall efficiency for convering rest mass into radiation, and t^ = Mc /L^ = 
4.4 X lO^yr is the Eddington time. 

However, even in the presence of a stabilizing external potential, an initially cold gaseous 
disk is prone to numerous instabilities (Binney & Tremaine 1987). Local stability to ax- 
isymmetric perturbations requires that the Toomre ratio Q = CsK,I{'kGYi) > 1, where 
E ~ (M/^ttR ) is the disk surface density [( ~ 0.1 is the ratio of disk to halo mass) 
and K ~ 2[GM/ R ) ' is the epicyclic frequency. For the present system we have Q ~ 
3(^/0.1) Mg i?3 (cs/10kms~ ), indicating marginal stability. Depending on the poorly 
known details of the thermal processes involved, some disks could easily become locally un- 
stable everywhere. In such a case we expect further fragmentation of the disk material into 
cold gas clouds or stars (as observed in the disks of present spiral galaxies), and the growth 
of a central supermassive black hole by viscous accretion would not take place. 

Even if local stability is achieved, a self-gravitating disk can still be unstable to global, 
non-axisymmetric modes. In the absence of a central compact mass, the self-gravity of the 
disk (whose mass interior to r scales like r close to the center) will always become dominant 
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over the external halo gravity (mass oc r ). In that case the central region of the disk would 
behave like an isolated, self-gravitating disk of cold gas, which is known to be violently 
unstable to non-axisymmetric modes (Binney & Tremaine 1987). To obtain global stability, 
we require the presence of a central seed mass Mg such that for all r. 



47r „ GMs 

3 -r.- . ^2 



GS < —Gphr + -^, (3) 



where the halo density ph ~ M / {Att R^ / 3) . With E ~ CMI{'kB?) the stabihty condition (3) 
translates to the constraint 

Ms > 10^ Mq X Mg (-^^ . (4) 

For ( ^ 0.1, the above inequality is consistent with the fact that the empirically deduced 
black hole masses are always ^ 10 M© (cf. point 4 in the introduction). 

The presence of a sufficiently massive central seed object is essential for the growth of a 
massive black hole near the center of the gaseous disk. One possibility is that the central seed 
forms during the initial collapse out of gas with low angular momentum (Eisenstein & Loeb 
1993). As illustrated in Figure 4, there are two routes that could lead to the formation of a 
seed black hole out of this gas. If the cooling time of the gas is shorter than its viscous time 
(cf. Eqs. [1] & [2]), the gas will settle into a self-gravitating supermassive disk, supported 
mainly by rotation (cf. §4.2). Angular momentum transport by nonaxisymmetric instabilities 
or by viscosity could eventually lead to the formation of a seed black hole for a sufficiently 
compact initial disk. On the other hand, if the viscous time is shorter than the cooling time 
the gaseous cloud could collapse quasi-spherically. The collapse would either lead directly 
to a black-hole (e.g., via radiative viscosity; cf. Loeb & Laor 1992) or could be slowed down 
by radiation pressure if the entropy is high enough. In this latter case, a supermassive 
star would form as an intermediate state (Wagoner 1969; Shapiro & Teukolsky 1983). The 
dominance of radiation pressure in both supermassive stars and compact supermassive disks 
would prevent fragmentation and star formation in these systems. 

4-2. Supermassive Stars and Disks 

It is possible that during the initial collapse a small amount of gas with very low specific 
angular momentum will condense near the center of the potential well and collapse quasi- 
spherically. As its density rises, the gas must eventually become optically thick. When the 
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photon diffusion time becomes longer than the gas infall time, radiative coohng becomes 
inefficient and the gas heats up adiabatically, with the pressure rising a.s p ex. p . If radiation 
pressure becomes dominant during the collapse, the adiabatic index of the gas will approach 
the value T = 4/3, and pressure support will never be able to halt the collapse since dp/dr oc 
1/r oc GMp/r . However, if the entropy production in shocks is sufficiently large, significant 
pressure support may build up, leading to the formation of a supermassive star (Hoyle & 
Fowler 1963) ^ 

In most cases, however, rotation is likely to affect the collapse earlier than the pressure 
support. Without loss of angular momentum the rotational velocity of a fluid element scales 
as f0 oc r~ , while its radial infall velocity rises only as r~ ' . Therefore a supermassive 
star of radius i?sm can form out of a cloud with an initial radius Rq only if the initial ratio 
of rotational to radial velocities is <C {Rsm/ Ro) ■ While this constraint may be satifled in 
rare cases (Eisenstein & Loeb 1993), a majority of systems must lose angular momentum 
very efficiently for the flnal stage of the collapse to be quasi-spherical. In the absence of an 
efficient mechanism for removing angular momentum, a supermassive disk (Wagoner 1969) 
must form. 

It is usually assumed that supermassive stars are convective, and can therefore be mod- 
eled as polytropes with an adiabatic index T ^ 4/3 (Wagoner 1969; Shapiro & Teukolsky 
1983). Since we could not flnd a justiflcation for this assumption in the literature, we pro- 
vide a simple but rigorous proof in appendix. Following equations {A2) and (^3), the star 
radiates at the Eddington luminosity, Le = 4:TrGMmp/a^ = 1.3 X 10^^ ergs"^ (M/IO^M©). 
Without rotation the star would contract on a timescale much shorter than the Eddington 
cooling time, t^ = Mc / Le = 4.4 x 10 yr, until it gets to a radius where the general rela- 
tivistic instability sets in (Shapiro & Teukolsky 1983). A rotating star must shed mass and 
angular momentum in order to reach this state. Because of the mass loss, and since rota- 
tion has a stabilizing effect, the radius of the star at the onset of gravitational collapse will 
be smaller than the value 1600(GM/c^) X {M/IO^MqY^'^ obtained in the nonrotating case 
(Shapiro & Teukolsky 1983). For a uniformly rotating supermassive star, the mass shedding 
limit is reached at a very small value of the ratio of rotational kinetic energy to gravitational 
binding energy, (r/|H^|)niax ~ 0.007 (Zel'dovich & Novikov 1971). In reahty, however, the 
evolution of a rotating supermassive star will depend sensitively on the internal distribution 



5 The question of entropy production in spherical collapse will be addressed in detail elsewhere (Thoul, 
Rasio & Loeb 1993). 
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of specific angular momentum. Even a small amount of differential rotation could stabilize 
the star against mass shedding, allowing much highers values of r/|H^| to be reached (see 
Bodenheimer & Ostriker 1973). 

Under the assumption that convection leads not only to uniform specific entropy, but 
also to uniform rotation , Bisnovati-Kogan, Zel'dovich & Novikov (1967) have calculated 
the evolution of a rotating supermassive star, modelled as an n = 3 polytrope. Once 
the contracting star reaches the mass-shedding limit, an equatorial outflow must appear. 
The subsequent evolution was modeled as a sequence of maximally rotating (r/|H^| = 
(T/|H^|)niax =constant) polytropic conflgurations with decreasing mass and angular mo- 
mentum. It was found that rotational support increases the lifetime of the supermassive 
star by up to a factor of 10 compared to the nonrotating case. Depending on the initial 
value of the angular momentum, between about 1% and 30% of the total mass is lost to an 
equatorial disk before the central star becomes unstable and collapses to a black hole. These 
results are for an object of mass 10 Mq, the only value considered by Bisnovatyi-Kogan et 
al. (1967). For a smaller star, it is possible that most of the mass will be lost to an outer 
disk. Therefore, even a very slowly rotating supermassive star could transform itself into a 
supermassive disk before becoming unstable to quasi-spherical collapse. 

Clearly, supermassive disks seem much more likely to form than supermassive stars. Un- 
fortunately, both the dynamical properties and the secular (viscous and thermal) evolution 
of isolated self-gravitating disks are not well understood theoretically. They depend sensi- 
tively on the unknown mass and angular momentum distributions, on the physical origin of 
the effective viscosity (which affects the heating rate and mass transfer rate as a function of 
radius and time), and on the transport of mass and angular momentum by dynamical insta- 
bilities. The secular evolution of isolated massive disks was discussed recently by Eisenstein 
& Loeb (1993), assuming a simple a-viscosity. The viscous time was found to be shorter 
than 10 yr (comparable to the star formation time) for the initial collapse of rare, low- 
angular-momentum objects in the universe. Most of the seed black holes were found to form 
just above the cosmological Jeans mass ~ 10 Mq (with roughly one seed object per bright 
galaxy), which, by coincidence, satisfles the constraint of equation (4). 



6 This assumption is apparently not satisfied in the convective envelope of the sun, where helioseimological 
data shows evidence for differential rotation (Brown et al. 1989). The solar rotation profile can be 
linked to the kinetics of convective turbulence (Kumar, Loeb & Narayan 1993). 
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5. SUMMARY 

We have performed numerical simulations of the hydrodynamic collapse of rotating clouds 
of cold gas in the centers of protogalaxies. We considered two simple limiting cases for 
incorporating angular momentum in the initial configuration. The first is solid body rotation, 
representing the well-organized fluid motions that could result from earlier dissipational 
evolution. The second is an irrotational velocity fleld, of the type produced by dissipationless 
tidal interactions with other distant masses. In both cases we flnd that most of the gas 
fragments into dense clumps. We model these clumps as a coUisionless system that interacts 
with the remaining gas only through gravity. During the evolution, a few subsystems form 
and then rapidly collide and merge together into a single virialized structure resembling 
the bulge of a spiral galaxy. The remaining gas forms a smooth thin disk that contains a 
fraction ^ ;^ 0.1 of the total mass and is supported by rotation in the external potential of 
the coUisionless component. A central seed black hole of mass ^ 10 Mq X Mg[(/0.1) is 
necessary in order to stabilize the central region of the gaseous disk, allowing the smooth 
growth of the central seed by accretion. This minimum seed mass is consistent with the 
lowest values of the empirically determined black-hole masses in active galactic nuclei. A 
seed black hole could form directly through the initial collapse of gas with very low speciflc 
angular momentum. In some cases, a quasi-spherical supermassive star may form as an 
intermediate state. More commonly, rotational support will precede pressure support and 
lead to the formation of a supermassive disk as an intermediate state. 
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paper was written. Computations were performed on the Cornell National Supercomputer 
Facility, a resource of the Center for Theory and Simulation in Science and Engineering at 
Cornell University, which receives major funding from the NSF and IBM Corporation, with 
additional support from New York State and members of its Corporate Research Institute. 
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Appendix 

In the following we show that the radiative cooling of a supermassive star inevitably leads 
to the development of convection. Consider a supermassive star which is purely radiative 
initially, i.e. has a stable entropy profile. Without a nuclear energy source (nuclear burning 
can be ignored for masses ^ 10 Mq] cf. Wagoner 1969), the entropy of each electron-proton 
fluid element must be changing according to the local radiative heat flux i^ at a rate given 

by 

rup ot p 

where p, Tg and s are the mass density, temperature and speciflc entropy, and nip is the 
proton mass. If the opacity is dominated by Thomson scattering, the local radiative heat 
flux F is related to the radiation pressure gradient by 

^=-^-Vp„ {A2) 

o-j, p 

where a^. is the Thomson cross-section. Ignoring gas pressure and rotation, the hydrostatic 
equilibrium equation is simply 

-Vp^ = I {AZ) 

P 

where the gravitational fleld g obeys Poisson's equation, 

V ■g = -^-KGp. (yl4) 

Combining equations (yl2)-(yl4), we see that the RHS of equation (Al) is constant, 

1^ ^ 'i-KGrrij, , , . 

-V • F = = constant. {A5) 

p o-j. 

By taking the gradient of equation (Al) we then obtain 

d -. 2TrGmlVTa 

The temperature gradient in the radial direction is negative since heat flows out of the star. 
Therefore, a negative radial entropy gradient must develop eventually, leading to convection. 
For a rotating supermassive star, the above result holds as long as the rotation period is 
much longer than {Gp)~ ' . 
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FIGURE CAPTIONS: 

Fig. 1: Dynamical evolution of a spherical, uniformly rotating cloud. Gas particles are 
shown on the right, coUisionless (high-density) particles on the left. The units are defined in 
§2. 

Fig. 2: Time evolution of the fractional mass in the smooth gaseous component during 
the dynamical collapse of a uniformly rotating spherical cloud (solid line) and an irrotational 
ellipsoidal cloud (dashed line). 

Fig. 3: Dynamical evolution of an ellipsoidal, irrotational cloud. Conventions are as in 
Fig. 1. 

Fig. 4: Possible evolutionary paths leading to the formation of a central seed black hole. 
As shown in Figs. 1 and 3, the initial collapse of a protogalactic gas cloud results in the 
formation of a thin disk supported by rotation in the external potential of a coUisionless, 
bulge-like component. The inner region of the disk is initially unstable and may become 
optically thick. If the viscous time for angular momentum transport is shorter than the 
cooling time, the gas will become radiation-pressure supported and could form a supermassive 
star. The star contracts quasi-spherically and spins up as it radiates at the Eddington rate 
Le- It becomes unstable to mass shedding when the ratio r/|H^| ^ 0.01. If most of its 
mass is not lost to an outer disk, the star will ultimately become gravitationally unstable 
and collapse to a black hole. If, initially, the viscous time is longer than the cooling time, a 
supermassive disk would form. Such a self-gravitating disk is unstable to nonaxisymmetric 
instabilities when r/|H^| Ri 0.3. Through various sources of effective viscosity, the central 
part of the disk will lose angular momentum and may evolve directly into a black hole on 
the viscous timescale. 
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